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1. Inbroducbion 


The program of study described in this report addresses 
two broad topics. Our first objectiv'e was to apply the most 
powerful spectral and radiative transfer diagnostic techni- 
ques to the S055 data on magnetic loops and dark absorbing 
structures, so as to place new limits on physical conditions 
in the plasma. Our second; aim was to investigate, through 
numerical modelling of magnetically guided acoustic shock 
waves , to what .extent slow mode shocks could explain the spicule 
and fibril phenomenon, and elucidate the injection of mass 
and energy into the solar corona. 

One paper, co-authored by P. Poukal and J. Raymond is 
entitled "On the Thermodynamics of Magnetic Loops in the 
Solar Corona. " It describes the reduction of S055 spectra 
of coronal loop structures, and their ‘nterpretation in terms 
of temperature, density and thermal equilibrium diagnostics 
in the loop plasma- The main thrust of this paper is to dis- 
cuss critically to what extent high radiative losses in the 
coronal plasma, are to be interpreted as evidence for enhanced 
magnetic heating of these volumes. 

A second study focusses on the detailed physical 
conditions in the remarkable dark absorbing clouds visible 
in the corona on S055 EUV raster pictures, but not generally 
visible in optical or CaK photoheliograms. This work 
utilizes the large non-LTE computer code developed at 
SAG by B. Avrett, to study radiative transfex* in the solar 
atmosphere. A paper is being prepared for submission with 


J. Raymond and U. Avj.'Ott at the Center for Astrophysics* 

J. Raymond has acted as a paid consultant to AER, Inc., on 
these first two studies. 

A third paper entitled "Propagation of Magnetically 
Guided Acoustic Shocks in the Solar Chromosphere" discusses 
our basic approach to handling acoustic shocks in a specific, 
inhomogeneous model of the chromosphere and magnetic geometry. 
The main result is that even in an inhomogeneous static 
chromosphere, it is difficult to transmit acoustic shocks 
with sufficient efficiency, to heat the corona. This paper 
has been accepted for publication in "Solar Physics. " It was 
written in collaboration with M. Smart, who was employed at 
AER, Inc., under this contract during the suromer of 1979, 
after completing his senior honours thesis at Harvard on this 
subject, under my direction. 

A fourth paper, entitled "A Dynamical Model of Chromo- 
spheric Spicules, and Their Influence on the Mass and Energy 
Balance of the Solar Corona" is in preparation. This paper 
extends the results of the Poukal and Smart study to a moving 
atmosphere, and develops a specific model of spicule and fibril 
formation. This work has been carried out with the assistance 

of A. Munshi, a Harvard student employed at AER, Inc. in a 

\ 

part-time capacity, under this contract. 

In section 2 the data reduction, modelling and results 
are described in more detail. In section 3 we draw some 
useful conclusions from the study as a whole, and suggest 
several fruitful directions for future work. 
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Data Uecluction» ^rhGorGUxcal MoclolX^intj and RgsuIUs 
a) Thormodynainics of Magnetic Loops in tho Solar Corona 
i) Introduction 

I’he observation that coronal radiative losses from 
magnetic active regions are typically enhanced by an order 
of magnitude compared to the quiet sun, is often cited as 
evidence that heating, that is,* dissipation of mechanical 
or electrical energy, is higher in these regions. 

In fact, this deduction is only reasonable if the 
higher losses can be proven to originate in a volume of 
higher plasma temperature relative to the lower tempera- 
ture Tg of the surrounding corona. In that case, the second 
law would rule out the conduction or advection of heat from 
0?2 to T^, to balance the radiative losses. This would then 
necessitate a locally increased deposition of energy derived 
from dissipation of e.g. electric currents (Shklovskii 1965 
Gold 1967, Rosner et al 1978) Alfven waves (Alfven 1947, V7entzel 
197 8, lonson 1978) or channeled acoustic waves (Osterbrock 
1961, Poukal and Smart 1980). 

In px-actice, it has proven difficult to demonstrate 
the r©qxi.irements of thermal equilibrium and Tj_ > Tg, except 
in flares, where lines of e.g. PeXXIV demonstrate T^^ K 10 
>> T2* In general, the real error bars on plasma tempera- 
tures derived in active region loops and background corona 
(e.g. Pecker et al 1954, Piror and Zirin 1963, Mason 1975 
Gabriel and Jox-dan 1975) exceed the quoted differences AT - 


3 



In this study, wo bring to boar upon **his problem, a detailed 
analysis of the rich EUV line spectrum of several bright 
loops observed with the S055 spectrometer during the Al’M 
missions. To the best of our knowledge, this analysis is the 
first instance where S055 spectrometer pointings and wavelength 
scans have been located in well-defined loops at the limb. 

We apply emission measure analysis, and also line intensity 
ratios to study a) the thermal equilibrium of loop plasma 
b) the temperature distribution of plasma in the loop rela- 
tive to outside corona c) the density and element abundances 
in the loops. 

ii) Data and Reduction 

Through examination of the full S055 data base of 
spectra taken between June 1973 and January 1974, we were 
able to find three sets of spectra in which the pointing 
was definitely within a loop off the limb and not contaminated 
by disc features. The pointing information for the 5x5 
arc sec spectrometer aperture was obtained from the ATM 
atlas, and also from large-area S055 raster pictures in 
ultravio.let radiations, in which the aperture could be 
located directly, to within 5x5 arc sec. 

Fig. 1 illustrates sample rasters of loops, in ultra- 
violet lines covering the range of plasma temperature between 
1.5 X 10^ K (OIV554) to 1.6 x 10® K (MgX 625). The loca- 
tions of the spectrometer aperture are shown by small dark 


circles. 
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A total of 9 full spectra of high quality, wore taken 
at these pointings, covering in each case the spectral range 

O 

280 - 1340 A®, with 1.6 A spectral resolution. Our first 
step in reduction of the spectra, was to identify the various 

ultraviolet lines of interest, including several weak emis- 

. . 5 

srons rn the temperature range T > 10 K, never before used 

as solar plasma diagnostics. 

Calibrated intensities were then derived, for the 

lines listed in Table I, using the updated calibration of 

Reeves et al (1977) . 

iii) Results 

The emission measure defined as; 

/1.122T 

EM(T) =1 N„ dh 

yo.89lT ® “ 

v?as calculated, assuming constant gas pressure, from the 
absolute intensity in each line. The curves of EM(T) 
averaged over the spectra in each loop, are plotted in 
Rig, 2 for each of the three loops. For comparison, a fourth 
curve gives the emission measure distribution derived by us 
from the averaged quiet corona line intensities published 
from the S055 spectra by Vernazza and Reeves (1978) . 

As a measure of the error in the emission measure curves, 
we have used an equilibrium-excitation and ionization code 
(Doyle and Raymond 1980) to calculate the predicted line in- 
tensities from this emission measure distribution, for compar- 
son (Fig. 3) with the observed values. It is seen that 
agreement is generally quite good — to within about a factor 
two. 
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Toblo I 
—2 —I —1') 

Intensities (ergs era s sr of Lines Measured 


Ion 

X 

Jan 7 

Jan 8 

Jan 5 

cn 

1335 A 

244 

1.3 

366 


904 

4 



CIIX 

1176 

25 


93 


977 

66 

4 

150 

Nil 

1085 

15 


13 

NIII 

991 

10 


10 

NIV 

922 

4.5 




764 

32 

3 

25 

NV 

1240 

65 

38 

93 

Oil 

702 

17 

2 

20 

OIV 

789 

106 

16 

81 


554 

112 

13 

63 

OV 

760 

42 

7 

30 


630 

372 

75 

330 

OVI 

1034 

1135 

614 

1546 

NelV 

468 

19 




542 

18 


«... 

NeV 

570 

33 

4 

17 


481 

31 


31 

NeVI 

560 

49 

11 

60 


434 

207 

50 

247 

NeVII 

895 

21 

20 

20 


465 

259 

104 

268 

NeVIII 

775 

242 

222 

238 

MgVI 

401 

<590 

— — 


Mg VI I 

434 

286 

177 

209 

MgVIII 

436 

276 

235 

238 

MglX 

706 

15 

24 

33 

MgX 

613 

482 

621 

659 

Silll 

1206 

124 

mm 

67 

SilX 

950 

7 

16 

5 

SiXI 

581 

23 

42 

66 

SiXII 

505 

531 

850 

1092 

SVI 

937 

42 

23 

38 

SXIV 

425 

165 

248 

170 

FeXVIII 

974 

5 

7 

9 

FeXX 

805 
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Dovialiions oil bho ratio plottod in Pig. 3 from unity are a 
measure of primarily/ abundance errors and departures from 
ionization and excitation equilibrium in the loop plasma. 

Two main points emerge from the relative shape of the 
three loop curves in Pig. 2/ as compared to the background 
emission measure curve; a) the loop curves do not turn over 
at higher temperature log T i 6.5/ as rapidly as the back- 
ground curve. This suggests that there may be relatively 
more hot plasma in these structures than in the quiet 
corona, b) there is also a lot more cool plasma at. T < 10^ K, 
well above the limb in these loops, than in the quiet corona. 

Point (a) shows that even sunspot-type loops, visible 
most prominently in lines such as OIV, QVl and NeVll show 
a small emission-measure excess at coronal temperatures. 

i'/ 

This certainly supports the view that loops in general 

6 

exhibit regions of relatively high losses from T > 10 K 
plasma radiations. However, we stress, and it can be 
seen from Pig. 2, that it is exceedingly difficult to 
demonstrate from emission measure analysis alone, that higher 
plasma temperatures are actually achieved in the loops than 
in the background. 

Table II shows the results obtained from the line inten- 

sity ratios, that we have used to indicate departure from 

thermal equilibrium in the plasma. The line intensity ratios 

4 5 

given span the , temperature range roughly 6x10 K to 7x10 K. 
Comparison of the values obtained in our loops, with the 
ratio calculated from our model for the case of thermal 
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equilibriunif gives an estimate of departure from thermal 
equilibrium. This difference is indicated by the disparity 
between the temperature at which the lines are actually 
formed/ and the temperature at which the ion should be 
formed/ if thejnnal equilibrium were rigorously maintained. 


TABLE 11 


Ion 

liine 

Ratio 

Loop 

Background 

Model 

^9 ^lines 

Log T . ^ 

^ xon 

NeVII 

1 (895) 
I (465 ) 

0.079 

0.045 

0.038 

5.40 

5.75 

NeVII 

1(562) 
1 (465) 

0.059 

- 

0,078 

<5.75 

5.75 

OIV 

I (554) 
1(789) 

1.05 

1.25 

1.24 

5.15 

5.20 

OIIl 

I (508) 
1(702) 

0.48 

0.43 

0.51 

4.90 

4.90 


In general, we see from the difference in the last tv/o 
columns , that the lines tend to be formed at lower tempera- 
tures, i.e. log log "^ion* consistent with 

the situation found in our model, when a plasma is allowed 
to cool rapidly from coronal temperatures. ^This supports 
the view (Foukal, 1976, 1978) that the cool loop plasma is 
coronal material condensing and cooling to chromospheric 
temperatures. This result is supported also by our prelim- 
inary analysis of several other line-intensity ratios , for 
which atomic parameters are less certain, and which we have 
thus not included in Table II. 
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b) Physical Conditions in Cool Coronal Clouds 


i) Introduction 

Optical observations in the and CaK lines show that 
large volumes within the hot (T > lO*’ K) coronal medium are 
occupied by much cooler T 10^ K plasma. This cool material 
is localized in discrete regions seen in and CaK as dark 
filaments on the disk, and as bright prominenc^iis on the 
limb. It has been shown in previous EUV studies (e.g. 

Schmahl et al 1974, Orrall and Schmahl 1976) that this 
material is optically thick in the Lyman continuum, at X < 

912 A. 

The hot (T > 10® K) coronal medium between these fila- 
ments and prominences is coh/entionally taken to be optically 
thin in the ultraviolet. This assumption is essential to 
physical analyses of solar and stellar coronae (e.g* Pottasch 
1963, Withbroe and Noyes, 1977) . 

But several analytics of ultraviolet data (Withbroe 1970, 
Kanno 1979, Foukal 1978, Schmahl and Orrall 1979) show evidence 
for important absorption by neutral hydrogen outside these 
discrete optical features. In a previous paper (Foukal 1978) 
we showed very dark, discrete absorbing structures in S055 
raster pictures taken in EUV lines and continue. In this 
study we have analysed additional SOS 5 ultraviolet data on 
these dark coronal "clouds," and also compared them with good 
quality optical data. We also used a non-LTE radiative 
transfer model to study the physical conditions in these 
structures, and relate them to coronal dynamics. 
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ii) paUa and Reduction 

data used in this study- war® polychromatic EUV 
raster pictures from the S055 spectrometer. The S055 rasters 
taken between Juno 1973-January 1974 were searched for 
distinct dark structures in EUV radiations. Pour examples 
of these dark clouds are given in Pigs, 4, 5/ 6, 7. We see 
distinct dark structures ranging morphologically between an 
arcade of parallel fibrils or low magnetic loops (Pig. 4) , 
an irregular “cloud” over plage (Pig. 5), an active region 
filament-like structure (Pig, 6 ) , and a time series of a 
loop-like structure (Pig. 7) . To examine the optical ap- 
pearance of these EUV dark clouds, we obtained the best 
guality photopeliograms and caK speotroheliograms from 
Big Bear and Sacramento Peak. These were scaled and oriented 
similarly to the EUV pictures, and are also shown in Pig. 4-6. 

We note that in general, there is no evidence for promi- 
nent dark optical structures, such as active region filaments, 
at the location of the distinct EUV clouds. Conversely, 
dark active region filaments (e.g. Pig. 5) seen prominently 
in are not always visible in the EUV radiations, or in 
CaK, although they can appear dark in some cases not illus- 
trated here. 

As v/e pointed out in our previous paper (Poukal 1978) , 
the opacity of EUV clouds at wavelengths such as X977, XI 03 2, 
X1335 located longward of the Ljrman continuum edge cannot be 
explained in terms of neutral hydrogen. The opacity at 
X1335 is likely to be caused by resonance scattering in the 
core of the strong resonance line of CII at that wavelength. 


10 


e 


However, at X977, X1032 wo require a eoureo oC eonUinuum opa- 
city never before considered in connection with the optical 
depth of the hot corona. This opacity source is CX, parti- 
cularly its 2p p continuum, with an edge at Xoo « llOOA. 

The rasters also yield useful information on the 
dimensions, height over the photosphere, and lifetimes of 
the EOV clouds. The basic observational constraints are 
summarized in Table XXI. Our objective in modelling these 
structures was to determine the simplest model that would 
CKplain these constraints in a physically reasonable and 
self-consistent manner, Xt is beyond the scope of our 
analysis to establish the uniqueness of this model, although 
we do compare our results with some of the closely studied 
models of active region filaments and prominences. 

In our modelling of those optically thick clouds, we 
used the non-LTE radiative transfer code "pandora" developed 
ac SAO for studies of the solar photosphere and chromosphere 
(Vernazza et al 1973). Our basic physical model is illus- 
trated in Pig. 8a. It is a plane parallel infinite slab 
lying in the xy plane, suspended magnetically above the 
photosphere at a height of 20,000 km, with the Z-axis normal 
to the local solar photospheric surface. The density p in 
the slab is taken constant, with magnetic stresses assumed 
to balance any pressure gradients arising from temperature 
gradients in the Z-direction. We then use the Pandora code 
to calculate the ionization structure, opacity in HI and 
Cl, can lines and continue, and emergent fluxes in various 
strong lines, such as Lya and CXI 1335, Ha and CaK. 
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Table III. 


Summary of Observational Constraints 


1 . Morphology in EUV radiations ; 


2. Characteristic Horizontal 
Scale; 


3. Characteristic Vertical 
Scale ; 


4. Characteristic Thickness in 
Vertica3, Direction: 


5. Optical Depth t(X): 


6. Brightness Contrast C(X) 
relative to surroundings: 


Distinct dark loops, 
filamentary structures, 
or irregular shapes. 


b ^ 5 X 10^ km 


h<l X 10^ km 


< 4 

A h 1 X 10 km 


T(X) 1 at; XX1335, 1216, 
1032, 977, 896, 625, 554. 

't(X) <1 at; Call K-line, Ha. 


C(X) c< 1 at XX 1335, 1216, 
1032, 977, 896 and 554 
C'(X) '^1 at X625, Call 
K-line, Ha. 


iii) Kosults 

After several adjustments of the model primarily dictated 
by the need to produce high column density of Cl, but low 
column densities in Call and in the n«2 population of HI, we 
have found a temperature profile T(Z), and density in the 
I slab, illustrated in Pig. 8b. The temperatures gradient at 

* A 

T > 3 X 10^ K is very large, and is determined by the constraint 
that transition region emissions (e.g. Clll 977, OlV 544, OVl 

I 

1 1032) are very weak, even though the density in the very 

& 

\ dark structures is very high (Dupree, Poukal and Jordan 1976) . 

i 3 

Overall, our calculations show a very cool (T 6 x 10 

K) gas in the central regions, with an exceptionally steep 

r temperature gradient to coronal temperatures. Comparison 

I with conditions in prominences (e.g. Orrall and Schmahl 1976) 

and active region disc filaments (e.g. Tandberg-Hansen 1971) , 

indicates that the EUV clouds are even denser than these 

optical structures, and the temperature jump from very cool 

plasma (T ^ 8 x 10^ K) to coronal values (T > 10*^ K) is 

much (at least a factor 10 at 1 x 10^ K < T < 3 x 10^ K) 

steeper. 

We suggest, from the comparable lifetimes of the EUV 
clouds, and their diverse morphology, that they are not a 
particular and separate class of coronal structures. 

Rather, we believe them to be simply situations in which 
the isobaric condensation (e.g. Kleczek 1958) of coronal 
matter in filaments, fibrils, and loops has achieved excep- 
tionally high densities and corrospondingly vbry low tempera- 
tures. 
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c) Propagation o£ Magnotieaily GuicJod Acoustic Shooks 
in the Solar Chronio sphere 

i) Introduction 

The roost attractive roechanisro for heating of the outer 
solar atwosphero above the toroporature roinimuro has been 
dissipation of acoustic shock wave energy. This roechanisro 
(Biermann 1946, Schatzroann 194 9) utilizes the observed 
turbulence at the solar photosphere to generate sound waves 
that then steepen to shocks. Both the generation of the 
sound and the propagation and dissipation of the shocks can 
be treated using relatively well-understood fluid properties. 

However, recent observations from the OS 0-8 satellite 
(Athay and White 1978, Bruner 1978) have shown that there 
is insufficient power (<^ 1 x 10* ergs cro ^ sec ’^) in .travel- 
ling waves at periods P ^ 40 secs, to balance measured coronal 
losses. This makes it unlikely that these long-period waves 
are responsible for coronal heating. At the same time, 
acoustic shock calculations (Ulmschneider 1971, Bird 1964b) 
have shown that in a plane parallel static atmosphere, short 
period (P 50 secs) waves will be dissipated in the chromo- 
sphere, well before even reaching the corona. Thus the 
OSO-8 observations at P K 40 secs, and this theoretical 
result for P 50 sec waves would, when taken together ap- 
pear to rule out coronal heating by acoustic shock waves 
of any period. 

However, observations sensitive to short acoustic 
wave frequencies P < 40 secs (Deubner 1976, Athay and White 
1979) do show possibly very large flux densities at these 
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hiyh frequencies, around the temperature ininimum. It then 
becomes of considerable interest whether the theoretical 
result derived for plane parallel, static chromospheres also 
holds for a more realistic (ie inhomogeneous, moving) chromo- 
sphere . 

In this study, we have integrated the shock propagation 
equations through a static, but inhomogeneous chromosphere. 
Specifically, we determine the propagation and dissipation 
properties of acoustic (slow mode) shocks guided along the 
field lines of the intense and rapidly diverging magnetic 
field in the chromospheric network. 

ii) Equations and Model Atmosphere 

The propagation of an acoustic shock through a gravi- 
tational atmosphere has been treated by the method of character- 
istics for several cases of relevance to the solar atmosphere 
by Bird (1964a) and by Kopp (1968) . A total differential 
equation is derived, relating the shock strength a - M^-l 
or Mach number M , to variations in the density p, temperature 
T, and flow cross-section A, ahead of the shock. 

For the special case of a periodic train of weak shocks 
in a static atmosphere, the shock propagation equation takes 
the form; 

1 ^ _ 1 r 1 dao . 1£ 1 dA 4a -v 

a dX 2 “ a dX a 2 “ A dX “ a P ^ 
o o o 

Here a, g, P, y represent respectively the sound speed, 
acceleration of gravity, shock period, and the ratio C /„ 

Ir 
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of spocifie heats. The subscripts o an<3 1 denote values 

ahead of, and behind the shock. 

A detailed consideration of the various terms on the 

right hand side has been given in our paper (Poukal and 

Smart 1980) . The basic procedure in this analysis was to 

calculate an initial shock strength cx^ at an initial height 

X s= Xq in the atmosphere, and then to integrate this equation 

numerically, using a finite difference scheme, taking a^CX) 

as given by a chromospheric model e.g. (Vernazza, Avrett 

and liOeser 1973, Kopp 1968) and using a model of the magnetic 

field geometry, to determine 1 dA . 

A dX 

Our model of the magnetic field lines was obtained by 
solving the magneto-static equilibrium relation; 

Vp = J X B , 

// 

!/ 

,/ 

,?which was reduced to a balance between the external gas 

( 2 

Ipressure p(h) and the magnetic pressure B (h)/87r, at every 

hV^ight h above X^. The profile of A(h) is given in 

Table IV. 

The calculations were performed both at constant y » 

5/3, and with y variable. The atmospheric model of T, p 
and y that we used in our integrations is given in Pig. 9, 
iii) Results 

To check the treatment of shock wave propagation and 

dissipation used here, we first compared our results for a 

plane parallel (1 ^ - 0) /^atmosphere with Ulmschneider ' s 

A dX iJ 

(1971) previous calculations. For these runs, we used values 
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TABLE IV 


Magnetic Geometry 


Height (Jem) 

f 

Radius (Jem) 

2 

Area (Jem 

below 750 

170 

9.4x10^ 

at 900 

690 

1.5x10*" 

1200 

1730 

9.4x10*" 

■ < 1500 

2770 

2.4x10*^ 

1800 

3810 

4 . 6x10 

2000 

4500 

6. 4x10”^ 

above 2000 

4500 

6 .4x10 
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of and similar to his computations, in which the 

shock propagation and dissipation were treated quite differently. 
The results are given in Pig. 10, for shocks of 15 and 50 
secs period. It can be seen that the agreement is good, 
namely, the flux density remaining in the shock as a func- 
tion of height is closely similar. If anything, Ulmschneider ' s 
code predicts higher losses than ours. 

In Fig. 11, we show 'che losses suffered in a wave 

II I 

propagating up a div^irging magnetic flux tube. Here the 
relative losses are much less. Even the 15 sec wave, the 
most heavily attenuated, retains fully 70% of its flux 
density to 1500 km above the height of shock formation. 

However the results of Fig. 11 are obtained at relatively 
low initial shock strengths a^. Although the transmission 
is high, the total fluk, integrated over all flux tubes in 
the network, is still too low by a factor 70, to heat the 
corona. The reason is that the fractional area of the photo- 
sphere covered by the flux tubes is only (e.g. Stenflo 
1976) about 1 x lO"^. 

Pig. 12 shows what happens if we attempt to make up 
for the small fractional area,’ by increasing the acoustic 
flux density (ie a ) in each magnetic tube. Once again, 
the transmission drops, so that the transmitted total 

- - 

flux emitted globally into the corona, is always too small / 

I 

to balance its observed losses. 

Our conclusions are; a) diverging magnetic fields may 
well redistribute acoustic flux dissipation vertically. 
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rolaUivc to non-magnetie xJhotospluQro and this could bo a 
factor in determining the vertical temperature profile of 
network magnetic flux tubes, and also plages, b) the global 
acoustic flux that penetrates even an inhomogeneous , static 
atmosphere is unlikely to be sufficient for heating the 
corona. 

d) A Dynamical Model of Chromospheric Spicules, and of 
Their Influence on the Mass and bnergy Balance of the Solar 

Corona. 

i) Introduction 

The most promising mechanism for accelerating spicules 

was originally put forward by Osterbrock (1961) and Uchida 

(1961) . It associates spicules and active region fibrils 

with the perturbation in temperature, density and velocity 

behind a slow-mode (acoustic) shock guided through the 

chromosphere by magnetic fields of e.g. the network. Wentzel 

and Solinger (1967) showed that in a uniform vertical magnetic 

field, a shock of plausible initial strength can produce 

temperatures and densities close to spicular values, when 

the gas velocity behind the shock is roughly 20 km sec , 

as observed in spicules. But that study did not consider 

the important influence of rapidly diverging magnetic field 

lines on the shock propagation. More generally, no attempt 

was made to explain the most basic and troublesome observation, 

namely; how does such cool and dense gas get so high into 

the corona, well above the ballistic trajectory h calculated 

2 

from the observed velocities, that is, h >> u • 




In this study, we use the shock propagation equations 
described in sect (c) above, to construct a reasonable dynami- 
cal model of the spicule phenomenon. 

ii) Equations, Physical jModel, and Results 
To calculate the propagation of the acoustic shock 
along diverging magnetic field lines, we used a similar propa- 
gation equation to that described in sect (c) above, but 
without the dissipative term, of the form 4a/aQp, that 
describes damping of the shock by its own wake. The reason 
for neglecting this term is that its ratio to other terms in 
the equation scales as 1/pf Where P is the shock wave period. 
For the high frequency shocks discussed in (c) , the term is 
important. But for long period disturbances (P ^ 300 secs) 
that would generate spicules, the term is negligible. As a 
first step in the actual calculation, we derive an initial 
Mach number at the height h ^ 750 km, where Ulmschneider ' s 
(1971) simple wave computation indicate the shocks will form. 
Judging from the velocities observed in granulation, this 
Mach number will be in the range 1.1 - 1.8. 

The propagation equation is then integrated upward 
through the diverging magnetic field geometry that was calcula- 
ted from magneto-hydrostatic equilibrium in sect (c) above 

(see Table IV). Using the values of M„ (h) from the propaga- 

s 

tion equation, it is then a simple matter to calculate the 
gas flow velocity u^, sound speed a^^, and density 
immediately behing the shock, using the standard jump 
relations (e.g. Landau and Lifshitz 1959). 
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The results are given in Pigs, 13, 14, 15. We note 
in particular, that the shock Mach number M_ drops around 
h 1000 km, due to the diverging field lines, and then re- 
covers when the field lines become parallel. The gas 

immediately behind the shock is accelerated powerfully, to 

-1 4 * 

outward velocities that reach 65 km sec at h « 1 x 10 km, 

behind the stronger of the two shocks shown. The gas tempera- 

5 

ture, calculated from the sound speed a^^, is about 1,6 x 10 K» 
and the particle density is 1.2 x 10 cm *"* 

We then calculate the cooling time of this plasma, 
after the shock has passed. The energy equation is; 


'“P H ‘i'W 

Here ^(T) is the radiative loss function for plasma of solar 
elemental abundances, as given e.g. by MeWhirter et al (1975). 
This yields for the cooling time, t , from the temperature 
T, behind the shock, to the observed spicular temperature 



If we assume, reasonably enough, that the plasma will cool 
isobarically , we can evaluate this integrali^to yield; 



9.5 X 10 


-10 


r 

B 


TdT "V 10 secs 
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At UhG sarao tiniGr the dccoXoration time under solar 
gravity, of this cooling gas is? 

tp “ Uj^/g 250 secs 

We thus have shown that the gas will cool rapidly to 

a temperature v/here it will become visible in and other 

optical radiations. The decelerated, cooled and condensed 

gas that is seen in these radiations as a spicule will have 

^*1 

an outward velocity of the order of 10 km sec""^, a tempera- 
ture T 'V' 1.5 X lO"^ and a density cm“^, as 

s 

observed (e.g. Beckers 1972). 


3« Over all Conclusions and Sucrgostions for Fugbhog Wogk 


The data analysis and theoretical modelling performed 
under this contract have led to a number of interesting con- 
clusions, that advance our knowledge of coronal dynamics. 

For one, our calculations of acoustic shock guiding 
have demonstrated how diverging magnetic fields in the net- 
work and plages can significantly increase the local trans- 
mission efficiency of acoustic power within the intense flux 
tubes. This opens up a new physical mechanism for heating 
chromospheric plages, that uses direct observations of vel- 
ocity fields in the photosphere, direct observations (e.g. 
Giovanelli, 1980 private communication) of field line diver- 
gence, and relatively well understood properties of acoustic 
shocks. On the other hand, our calculations also showed 
that oven in an inhomogeneous model of the chromosphere, it 
appears unlikely that the global transmission efficiency of 
high frequency acoustic power is sufficient to heat the corona. 
In this sense, our calculation strengthens the result (Athay 
and White, 1978) that acoustic shocks are not the primary 
heating source of the corona. 

Another clear result of our shock wave investigation, 
is that an exceedingly attractive model of spicule tempera- 
ture, density and velocity, can be calculated, using the 
acceleration of low coronal material by an acoustic shock of 
plausible (weak) strength, and then comparing its cooling 
time to its deceleration time under gravity. 
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Our study of shock waves guided along magnetic field 
linos suggests a number of fruitful directions to follow in 
the future. For one, it would bo a relatively simple matter 
to use our shock wave dissipation functions to calculate the 
relative temperature profile above the temperature minimum 
in the chromosphere, in a region of diverging fields, and 
in a non-'magnetic region. It would be interesting to compare 
such a profile with empirical models of faculae and plages 
(e.g. Chapman, 1970, Shine and Linsky, 1974), 

Another question posed by our results is what are the 
propagation characteristics of high freguoncy acoustic shocks 
in the moving medium represented by upwardly accoloratod 
spxcuies? The framework for such a calculation xs available 
in the propagation equations we have used. However a treat- 
ment including motions ahead of the shock requires explicit 
evaluation of complex non-adiabatic terms. We hope to address 
this question in the future, using accurate cooling rates for 
chromospheric gas, presently being calculated by Avrett at 
SAO. 

The work wo have performed under this contract on plasma 
temperature and density in loops has also laid the basis for 
a critical study of the thermodynamics of these loops. We 
have completed analysis of all the S055 spectra relevant to 
the investigation, carried out the emission measure analysis, 
and also most of the analysis of line intensity ratios . Our 
results indicate that there is indeed evidence to believe 
that enhanced energy deposition takes place in loops, since 
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the plasme temper at uire does seem to achieve higher values 
than in the background. We now are running final model cal- 
culations/ using now atomic croa* -sections rocently made 
availablo from the Belfast group, to comploto tho lino inton- 
sity-ratio analysis, before publishing our oonclusions. Wo 
expect that our results wrll yield the first proper demon- 
stration that active regions are in fact "hottor", and point 
out tho thormodynsimic consequencos of this result. 

Finally, our analysis of cool absorbing clouds has 
H lelded a physically consistent picture of the density and 
temperature conditions in those, remarkable structures. Our 
identification of these structures in the S055 rasters hau 
given now evidence for tho existence of largo opacity by 
cool gas in the coronal medium, established the basic physical 
conditions in these structures, and laid the basis for even- 
tual understanding of their dynamics. The recognition that 

O 

opacity is present even longward of 1912 A, has important 
implications for the use of plasma diagnostics in loops and 
other active region coronal structures. The close proximity 
of the ultra-cool material to sunspots confirms o\ar previous 
evidence that the corona over sunspots is in general a region 
of anomalously low temperature, so that the coolness of an 
umbra extends right through the photosphere, chromosphere 
and to at least 50,000 km in the corona. 

We acknowledge valuable assistance in our work from E. 
Avrett at SAO, whose calculations using the Pandora code have 
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been of great assistance in analyzing conditions in the cool 
clouds. Also R. Moore at BBSO and S. Hinata at SPO have 
helped us in providing photoheliograms and CaK spectro- 
heliograms that wq have used in that same analysis. 
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ultraviolet raster pictures of loops 
observed with the SOS 5 spectrometer. The 
format of these mini-rasters, taken to pro- 
vide pointing information for the spectra, 
is about 2 ' X 5 ' . In the vertical column 
on the left (January 5th) , a loop structure 
and pl\ame are seen in McMath 12690 on the 
SE limb at GMT 1711, in om554, OVIM032 
and Mg XX625. The two top pictures in the 
right hand column show a bright plume over 
a spot in McMath 12694, at the E limb on 
January 7th GMT 14:20, in NeVIlX465 and 
MgVIIIX436. A loop in the same active region 
on January 8th at GMT 01:13, is seen in the 
lower two frames, in the same two lines. 

Plot of emission measure against tempera 
ture for the three loops observed, and for 
"background" corona, over the temperature 
range 1 x 10*^ K < T < 5 x 10^ K. The lines 
used are listed in Table I, 

Plot of observed intensities in the 
lines used to calculate emission measures 
in Pig. 2, divided by predicted intensities 
of those; same lines, as calculated from the 
resulting plasma temperature distribution. 


4 Dark EUV absorbing niatorial near active 

region McMath 12503 on August 28, 1973. The 
location of the dark structuace is indicated 
by an arrow in the Lyman-alpha raster picture, 
taken at GMT 15; 01. The same feature is seen 
also on rasters taken simultaneously in Dyman 
continuum X896, and lines of CIIX1335, 
CIIIX977, OIVX554 and OVIX1032. The feature 
is not seen in MgXX625. Arrows also indicate 
the location of the EUV dark feature, on 
optical photoheliograms in (Big Bear photo- 
graph taken at GMT 15:13) and in CaK (taken 
at McMath-Hulbert at GMT 14:30). The optical 
pictures aro printed to similar scale and 
orientation as the EUV rasters. The feature 
is not seen in these optical radiations. 

Fig. 5 An irregularly shaped dark EUV absorbing 

structure in active region McMath 12507 on 
September 2, 1973. An arrow indicates the 
dark structure on the Lyman-alpha raster. 

It is also seen prominently on raster pic- 
tures, taken simultaneously at GMT 02:05, in 
resonance lines of CII, GUI, OIV and OVI. 

It is less visible in MgX radiation. No 
Lyman continuum observations were taken at 
this time« The closest available photo- 
heliograms were taken at Big Bear 2h 22 m 
before the EUV rasters, and 13h 21 m after- 
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Pig. 6 


I; 

\ 




f 


Pig. 7 


Pig. 8 


ward. The location of the EUV structure 
is indicated by an arrow in the earlier H 
photo. The closest available CaK was taken 
at Sacramento Peak llh 39 m after the EUV 
rasters. An arrow indicates the location 
of the EUV dark structure on this CaK frame. 

A dark EUV-absorbing structure in active 
region McMath 12535 on September 24, 1973, 

GMT 00:32. As in Pigs. 1, 2, arrows indicate 
the location of the structure in the Lyman 
raster picture and show the location of the 
same area also in and CaK photoheliograms. 
It is hardly visible in the MgX EUV raster. 

The two H pictures were taken lOh 02 m before 

v» 

(on left) and 8h 16 m after the EUV rasters. 

The CaK picture was taken at Sacramento Peak 

13h 27 m after the EUV rasters. 

Time sequence of a dark loop-like EUV 

structure near active region McMath 12474 on 

August 7, 1973. Arrow indicates the structure 

in the CIIM335 resonance line, at t^^ = GMT 

04:08. It is shown at t^ also in radiations 

o 

of OVI and OIV. Later frames of the sequence 
show the structure in these radiations 22m 
and 39m later. 

In (a) we illustrate the basic geometry 
of the sliib model. In (b) we give the "best 
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fit" temperature distribution T(Z) derived 
• from our iterations. 

Pig. 9 The vertical profiles of temperature T, density 
p and the ratio of specific heats Y/ the 

principal model atmosphere used in these cal- 
culations . The density profiles wore calculated' 
on the assumption of hydrostatic equilibrimn. 

In one case y was held constant, in the other 
case it was allowed to vary as shown above . 

Pig. 10 Curves of flux density P (h) remaining at a height 

h above Tq g =1, in a plane-parallel atmosphere 
of Y = 5/3. The curves shown here represent 
shocks of 50 second and 15 second periods. The 
curves designated U are from Ulmschneider (1971) . 

Pig. 11 Curves of flux A(h)P(h). remaining in a diverging 

• magnetic flux tube, for shock periods of 15 s and 

50 s. The initial flux values correspond to the 

same initial flux densities as ip Pig. 2, but are 

multiplied by the cross-sectional area 9.4x10^ 

2 

]cm of a network flux tube at h = 750 1cm. 




Fig. 12 curves of flux A(h)F(h) remaining in a diverging 
magnetic flux tube for a higher initial shock 
strength Uo “ 0.3. Shock periods of P = 15, 28 
and 50 s are shown. The curves given are calcu- 
lated with Y ® 5/3, thus for highest transmission 
The dashed curve was calculated with the model of 
Gabriel (1976) . 

Fig, 13 Plots of the shock Mach number against height, 
h, above t_ c “ 1/ for initial Mach numbers at 
h = 750 km of 1.4, 1.8. 

Fig. 14 Plots of the density immediately behind the 

shock, against height, h, for initial shock 
Mach numbers of 1.4, 1.8. 

Pig. 15 Plots of the sound speed a^^, and the gas velocity 

u^, directly behind the shock, against height, 
h, for two initial Mach numbers at h = 750 km; 


1 . 4 and 1.8. 
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